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This is the second in the series of papers on transport phenomena along random rough surfaces. 
We apply our simple general approach [lj to transport in very narrow channels, when the particles 
wavelength is comparable to the width of the channels and the motion across the channel is charac- 
terized by discrete quantum states. The discrete nature of the spectrum leads to a non-analyticity of 
transport coefficients as a function of thickness or channel width, especially for degenerate fermions. 
Surface inhomogeneity leads to both on-level scattering and interlevel transitions. As in jjj], trans- 
o\ ' port coefficients are expressed explicitly via the correlation function of surface inhomogeneities. The 

shape of the curves for the dependence of transport coefficients on the number of particles and/or 
film thickness is determined by the correlation radius of surface inhomogeneities and is not sensitive 
to their amplitude. For short range correlations, the curves assume a saw-like shape as a result of 
the interlevel transitions. With increasing correlation radius, the interlevel transitions loose their 
importance, and the saw teeth gradually decrease, reducing, in the end, to simple kinks on practi- 
cally monotonic curves. Applications include Gaussian and short-range 5-type correlations. Careful 
analysis of the transition from quantum to semi-classical and classical regimes allowed us to improve 
the accuracy of our previous classical calculations. 

PACS numbers: 72.10.Bg, 73.20.Fz, 73.50.Bk, 79.20.Rf 

m 

O ■ I. INTRODUCTION 

iy) ' Repeated collisions of (ballistic) particles with rough walls with random inhomogeneities result in gradual chao- 
^\ , tization of motion of particles along the walls, and lead to formation of the mean free path. Though this is a very 
old problem, important for different branches of physics, a simple consistent description has been absent. Recently 
we suggested a new simple quantitative approach to transport phenomena in systems with rough walls 0] (referred 
below as I; see also Q). Within this approach, one can easily express transport and localization parameters (such 
as mobility, diffusion, mean free path, localization length, etc.) for the motion of ballistic particles along the walls 
directly via the correlation function of wall inhomogeneities. The only restriction is that the inhomogeneities should 
be relatively smooth with the amplitude £ smaller than their correlation radius R and the thickness of the film (width 
q , of the channel) L, £ <C L, R . 

Our approach is based on the use of an explicit canonical coordinate transformation (similar to the Migdal transfor- 
mation in nuclear physics) which makes the initially rough boundaries flat. The price is a considerable complication 
of the bulk transport equation which, however, can be treated perturbatively. This transformation corresponds to 
stretching of the film (channel) width with the stretching parameters changing along the film in accordance with the 
exact shape of the boundaries. Since, in contrast to earlier, more heuristic approaches we are using an explicit 
expression for the coordinate transformation with the parameters given by the exact profile of the random boundaries, 
our reformulation of the transport problem with random rough walls as the transport problem with flat walls and 
randomly distorted bulk, is exact. The latter problem allows a straightforward perturbative solution which involves 
averaging over initial boundary inhomogeneities. 

In many cases, all the calculations after the coordinate transformation are really very simple and straightforward. 
In I we illustrated our method for a wide range of the problems. In this paper we will concentrate on the case of 
very thin films and narrow channels when the motion of (longwave) particles across the channel is quantized with a 
noticeable separation between the states. 

One of the effects that we want to describe is the following. In thin films with discrete levels for motion across the 
film, the increase in particle density and/or film thickness cause redistribution of particles between levels. In Fermi 
systems at T — ► this is a non-analytical step-like process. In this case, the lower levels are filled, while the higher- 
levels are unoccupied. The (gradual) change in the total number of particles results, at certain critical values, in the 
change of the number of occupied levels by one. These abrupt changes in the number of occupied levels should lead to 
a non-analytical dependence of the density of states and other thermodynamic functions on the number of particles. 
Obviously, this effect should also lead to singularities in the dependence of the transport coefficients on film thickness 
or density of particles. This saw- like effect has already been described in the case of bulk scattering la], and has 



1 



been qualitatively suggested in Ref. M for scattering by rough walls. We will be able to get quantitative results, and 
express explicitly the transport singularities via the correlation function of surface inhomogeneities. 

What is more, we will demonstrate that the picture is even more interesting than in [|||. We will see that the 
particle transport along the film is a non-trivial function of two parameters, NL 2 and L/R, where L is the thickness 
of the film with the 2D density of particles N, and R is the correlation radius of surface inhomogeneities. The saw- 
like dependence of the transport coefficients on NL 2 is well-pronounced only for short-range correlations of surface 
inhomogeneities, NR 2 <C 1. With an increase of the correlation radius R, the saw-like structure of the curves will 
gradually disappear, and the curves will become more and more smooth. It turns out that the saw-like shape of the 
curves is associated exclusively with the interlevel transitions the rate of which is determined by the parameter NR 2 ; 
if one artificially forbids the interlevel transitions even at small NR 2 , the curves will loose their saw- like character 
and will become similar to those at large NR 2 . 

Of course, the density dependence of transport coefficients becomes more and more smooth with increasing tem- 
perature even though the energy spectrum remains distinctly discrete. 

In the next Section we will present general equations for quantized transport in thin films with rough boundaries. 
Then, in Sec. 3, we will study transport singularities for degenerate fermions at T = 0. In Sec. 4 we will derive 
the corresponding transport equations for finite temperatures, and calculate the transport coefficients in Boltzmann 
temperature range. Our detailed analysis of transport of particles with discrete spectrum will demonstrate that the 
accuracy of previous calculations for continuous spectrum (i.e., in classical and/or semi-classical limit) in I can be 
improved considerably. The consistent transition from discrete to continuous expressions will provide us with much 
more accurate expressions for the squares of the ^-functions, which appear in the calculations of the squares of the 
matrix elements, than the approximations used in purely classical calculations I. The improved classical transport 
calculations will be described in the Appendix. 



II. QUANTIZED STATES ACROSS THE CHANNELS 

As in I, we will consider a film (channel) of the average thickness L with rough boundaries x — L/2 — z) and 
x = —L/2 + £2(2/) z ). The small boundary inhomogeneities, ^1,^2 L, are random functions of coordinates s = (y, z) 
along the boundaries, (£1) — (£2) = 0, with the correlation function 

Cife(|si-s 2 |) =te(s 1 )a(s 2 )), (1) 
Qk (q) = / d 2 s e^-/ R d k (s) 

where q is the momentum along the wall. In homogeneous systems, the correlation function depends only on the 
distance between points |si — S2 1 and not on coordinates themselves. In the absence of the bulk relaxation, the results, 
as in I, depend only on the function £ (s) = £1 (s) + £2 (s) and the correlation function £ |si — S2 1 = (£(si)£(s2)) = 
Cn + C22 + 2C12. 

Though we can calculate the transport coefficient for arbitrary correlation function £ (s), we will supplement general 
expressions by the most practical examples of Gaussian correlations of the surface inhomogeneities of an average height 
t, 

C (s) = £ 2 exp (s 2 /2R 2 ) , C (q) = 2irl 2 R 2 exp (-q 2 R 2 /2h 2 ) (2) 

including the limiting case of the very small correlation radius R, i.e. the <5-type correlations, 

( (s) = l 2 R 2 S (s) /s, C (q) = 2ir£ 2 R 2 (3) 

Note, that the condition I <C R does not mean that our approach is applicable to long-range correlations (large 
size inhomogeneities) exclusively. The scale for the effective correlation range in Eq. (|^) is defined by the particles 
wavelength A. For longwave particles A > one deals effectively with the short-range <5-type correlations (g), while 
in the opposite case of long-range correlations A <C R one should consider the full Gaussian expression (||) . 
Our approach is based on the use of canonical coordinate transformation 

L\x-\(t 2 (y,z)-Z, x (y,z))\ 
which corresponds to the following change of the form of the bulk Hamiltonian H = p 2 /2m: 
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(for details see I). The randomness of inhomogeneities, (£) = 0, leads to the randomness of the bulk "perturbation" V, 
(v\ = 0. Thus, the transformation (|j) reduces the transport problem between rough walls to an equivalent transport 

problem with ideal specular walls, i ff(L/2) — ^(—L/2) = 0, but with a distorted bulk Hamiltonian (Q). The latter 
problem can be treated in the same standard perturbative way as for any random bulk imperfections or impurities. 

In thin films, the motion of particles across the films is quantized, P Xj = nj /L, while the motion along the films 
can still be considered classically or semi-classically. Then the wave functions which correspond to the unperturbed 
Hamiltonian (H) with V = have the form 



*j (r) = y/2/v exp (iq • s) sin (njX/L) 



(6) 



(vq is the volume). In I we discussed the simplest situation when the distance between levels @ is so large that the 
transition between them are practically forbidden. Then the transport problem reduces to independent calculations 
of the mean free paths Cj for each level. The latter problem was treated as classical 2D diffusion on each level. 

The full quantum problem should include the transition between levels (^) with different j caused by the perturbation 
(|b|). Surface inhomogeneities lead to the change of quantized energy levels e^Q, interlevel transitions j — ► j' with 
diffusion between energy levels, and diffusion of particles on the same level. Since the perturbation is small and we 
are interested only in transport along the film, we can, as in I, neglect perturbation-induced corrections to the energy 
levels and neglect the distortion (||) everywhere except for the collision integral in the perturbative transport equation 
(i.e., neglect all distortions in dynamic, l.h.s. of the transport equation). Then the transport equation becomes a set 
of equations in the distribution functions rij coupled via collision integrals Lj : 



d t n (ej, q) H d T n (ej, q) + F ■ <9 q n (e/,q) = Lj {m} 
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(J) 



where tj - 
Wjj> (q, q')between states (j, q) — » (j',q')j 



(njh/L) + q 2 /2m. The collision integrals Lj are determined by the transition probabilities 
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The transition probabilities are given by the squares of the matrix elements of the perturbation (|^) : 

2tt 



33' (q,q') = y 



UfMMi)-4 (q')) 



(8) 



(9) 



Since all the equations include only the squares of the matrix elements of the "perturbation" V, the averaging over 
the random surface inhomogeneities does not cause any problems leading directly to the correlation function £ (s) . 
The matrix elements of the perturbation (|^) are 
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and the collision integral 

La = 



2mL 



L 2 



2mL 



■e(q-q') tf v -Qu) 



4J 



3 + 3 3 

for particles on each level j becomes, after averaging over inhomogeneities, equal to 
f 



2i:h 3 m 2 L 2 



J d 2 q' ( (q - q') ^ (rij, (q') - nj (q)) 5 {e jtq . 



J(q-q') 2 + (^ 




(11) 
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This collision integral includes both on-level scattering (diagonal terms with 8jf) and interlevel transitions (off- 
diagonal terms with 1 — Sjj> ) which are induced by surface roughness. 

In I we assumed that the interlevel transitions disappear (i.e., that the off-diagonal terms in ( pi] ) can be neglected) 
with an increase in the distance between levels. We solved the transport equation without interlevel transition with 
an additional assumption that the film is very thin and, therefore, the components of particle momenta across the 
film are larger than those along the film, (nhj/L) 2 ^> (q — q') 2 . The results below show that these assumptions 
are self-consistent for large correlation radii of surface inhomogeneities, NR 2 3> 1. For smaller correlation radii the 
interlevel transitions and momenta along the film cannot be neglected even in the case of large interlevel spacing. 
As we will see, the approximation I, though valid in important limiting cases, misses the most interesting features of 
transport through very narrow channels or thin films. 



III. SINGULARITIES IN TRANSPORT OF PARTICLES WITH DISCRETE QUANTUM STATES: LOW 

TEMPERATURES 

Changes in particle density and/or thickness of the film lead to redistribution of particles between discrete quantum 
states j. This redistribution between discrete states may lead to a non-analytic dependence of transport coefficients 
on particle density and thickness of the film. Of course, this non-analyticity is more pronounced for degenerate Fermi 
systems at T — > when continuous increase in the number of particles leads, at certain critical densities, to filling of 
new levels with higher and higher values of j. Below we will look at this non- analytic effect in some detail. 

At T = 0, the Fermi momenta of fermions for the motion along the film q^Pon each level j are given by the overall 
Fermi energy (chemical potential) ef as 

" = - ((l L 

while the 2D density of spin-1/2 particles on each level is 



a U)2 

Nj = ^-s (13) 

(for simplicity we assume that the effective masses of particles on all levels are the same). The chemical potential 
fi = ep should, as usually, be determined self-consistently by calculating the total density of particles N, 




Eqs. (fill|)-()l4|) in convenient dimensionless notations, 



L \ 2 2 2 

2me F I — , Zj = —NjL 2 , z = V Zj = —NL 2 , (15) 



can be rewritten as 



zj = v - j 2 , z = J2 z J (16) 

The number of occupied levels S for the given value of z (i.e., for the number of particles NL 2 ) is given by the integer 
part of v x l 2 (z), 

S (z) = Int yV\ (17) 



All the levels with the indices j > S are empty, Zj>s = 0. Summation of Eqs. (|16|) from 1 to S defines the number of 
occupied levels S and the dimensionless chemical potential v as functions of the number of particles z: 



S = Int [y/v\ = Int 



£ (S+1)(2S+1) 
S + 6 



(18) 



V{Z) =^ + \(S+l)(2S+l) 
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[For computational purposes, it is more convenient to start the calculations by defining the number of occupied levels 
S, and to determine the interval of the values of z and v, which corresponds to this number of levels, basing on the 
value of S]. The changes in number of occupied levels S = 1, 2, 3, 4, 5, 6, ... occur at z = 0, 3, 13, 34, 70, 125, ... 
(i.e., in the points z = S 3 - S (S + 1) (25 + 1) /6). 

At T = 0, we look for the solution of the transport equation ^ in the form 

FL 3 



n-j (q) 



,(0) 



cos 8i 



where Oj is the angle between the momentum q^ and the external force F. Then, after the integration of the collision 

the transport equation reduces to the 



integral ( |11| ) with the Gaussian correlation of surface inhomogeneities 
following set of S dimensionless linear equations in \j (qfj - 

1/2 r 2 



R 2 



S(z) 



2w 2 z,R 2 



L 2 



27T 2 Z,i? 2 



L' 2 



-zUFx -,4,- 



L 2 



+2 V (1 - %0 j 2 j' 2 exp \-t: 2 (y/zj - ^z~7) 2 R 2 /2L 2 



(19) 



/l 27T 2 ,/Z-Z-7i? 2 

XiM^ h' 1 ' T 2 



3 2Tr 2 ^z-z-TR< 
- xFi -, 2, ^ 




/I 2^ y /zjz?& 
2' ' T 2 



(we will not give here similar cumbersome equations for the correlation function of a general form C( c l))- The 
conductivity (mobility) of particles is given by the solution of this set of equations, 



(7 = (7 = 7 <J^j' = 

yy 22 / j uu 



and can be conveniently parametrized as 



r (i) 



r (j) 



e 2 L 2 ( L 
n*M 2 T' i? 



(20) 



(21) 



The functions $ (z) for four different values of R/L are plotted in Fig. 1 (R/L = 0.05) and Fig. 2 (R/L = 1; 3; 5). 
The singular points correspond to change in values of S from 1 to 2 to 3 to 4... at z = 3, 13, 34... 
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FIG. 1. Function $(x), Eq.(|2l|), for R/L = 0.05 (solid line); dashed line - the same function calculated without interlevel 
transitions 
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FIG. 2. Function Eq.@, for R/L = 1; 3; 5 
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Another way of parameterization of equations, similar to the one used in I, is based on definition 

27r 2 ^/zJzJR 2 A-K^fzJzJtNR 2 8tt 2 <fzjzp f R\ 2 
T 2 = z = z \x ) ' 

where the effective particle wavelength A 2 = 2ir/N. This equation redefines the function <f> (z, L/R) as $ (z, X/R) = 
$ (z, yfi\/2R). 

Dramatic difference in shapes of the curves in Figs. 1 and 2 for small and large values of R/L is explained by the 
role of interlevel transitions. If one simply neglects the inter level transitions (the off-diagonal terms) in the collision 
integral (|TT|) , then the set of transport equations ( |l9| ) will decouple into S independent equations. It is fairly obvious 
that in this approximation the conductivity (^0|) is almost always a monotonic function of z = 2NL 2 /ir, though the 
critical values of z, which correspond to the change in the number of occupied levels S, are still responsible for the 
singularities (kinks) in the curves. Therefore, the saw-like nature of the curves is caused by the interlevel transitions 
exclusively. 

For comparison, Fig. 1 (dashed line) and Fig. 3 give the function $ (z) for on-level scattering exclusively, i.e., when 
all the interlevel off-diagonal terms in the collision integral ( pd] ) are artificially neglected. Obviously, the curves with 
and without transitions always coincide as far as z < 3 and there is only one occupied level. The differences show up 
only at z > 3. 




0.1 

1 10 100 

FIG. 3. Function Eq.©, for R/L = 1; 5 caculated without interlevel transitions 

On the other hand, the importance of interlevel transitions is characterized by the parameter NR 2 ~ zR 2 /L 2 . 
Since zy — Zj = j 2 — j' 2 , the exponent 



exp -it 2 (y/zj - y/zy~) 2 R 2 /2L 2 = exp 



N 3 R 2 -JNyR 2 



always makes interlevel transitions to remote levels \j — j'\ ^S> L/R negligible. These exponents show that the 
interlevel transitions and the resulting mixing of adjacent levels are very important only for not very populated levels 
with 2ir 2 ZjR 2 / L 2 <C 1. Thus, the contribution of interlevel transitions is very noticeable only for relatively small 
values of R/L, and decreases exponentially with increasing R/L. For this reason, the saw-like character of particles 
mobility becomes less and less pronounced with increasing R/L. At R/L = 5 the saw nearly completely disappears, 
and there is practically no difference between the curves in Figs. 2 (exact calculation) and Fig. 3 (calculation with 
artificially neglected interlevel transitions). Note, that the curves calculated with and without transitions always 
coincide for a small number of particles z < 3 when only one level is occupied and the transitions are impossible for 
energy reasons. 

The same parameters, 2n 2 z 3 R 2 / L 2 = 4-kNjR 2 , determine the values of the hypergeometric functions in (n9J) . At 
x 2 < 1, 1F1 ((2n - 1) /2, n, -x 2 ) ~ 1, while in the opposite case x 2 > 1, i-Fi ((2n - 1) /2, n, -x 2 ) ~ (n - l)\J\/Trx n . 
Therefore, at large 2ir 2 ZjR 2 / L 2 ^> 1, one can not only neglect the interlevel transitions, but also the hypergeometric 
functions with n = 3 and n = 4 in the diagonal terms of the collision integrals (|l9|) in comparison with the one with 
n = 2. Only under these conditions one can justify a heuristic assumption made in I and recover the result (1.53): 
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3 = 1 

- 2 S(z) 



e 2 L 2 f R 

t^u 2 r' l 



(22) 



L 



,*4 



Lj AR 2 ^ lFl (l,2,-2ir 2 z,R 2 /L 2 ) 



For small correlation radii NR 2 <C 1, all the terms in ( |l9| ) are of the same order, the hypergeometric function 
1F1 {(2n - 1) /2,n,-2TT 2 Z: jR 2 /L 2 ) - iF x ((2n - 1) /2,n,0) = 1, and Eqs.(||) can be simplified as 



zj^ 
R 2 



\xi + <i--„r + \z) ) - 2 Xj f; (i - %0i 2 i' 2 



(23) 



Then 
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$f z ^-— y (24) 

' ' LJ 4i? 2 ^ (j4 + 3 Zj f/2 + 5z 2 /8 ) + 5 (S + 1) (25 + 1) /6 - j 2 



In the opposite limiting case when for all levels zjR / L 1, the intcrlevel transitions and higher-order hypergeo- 

'3 

, 2 ' 



metric functions can be neglected, iFi (1,2,— x 2 ) — > l/y'Trx 3 , and 



Note that Eq.(^5|) for large iVi? 2 can be improved near the critical values of z which correspond to changes in number 
of occupied levels S. With the appearance of a new level S, the number of particles on this level, zs, and, therefore, 
zsR 2 /L 2 are small even for large R/L, and the contribution of this level is zs/ \F\ (|, 2, — 2tt 2 z s R 2 /L 2 ^) ~ zs, and 

not ZS7T 1 / 2 (27T 2 z,sR 2 1 '£ 2 ) 3 ^ 2 as it is implied by Eq.(f25|). Only later, when the hypergeometric function will become 
small, 2-k 2 zgR 2 1 L 2 3> 1, the contribution of this highest level will recover the form indicated by (p5|). 

The argument of the exponents and hypergeometric functions can be also written as the ratio of the particle wave- 
length to the correlation radius of surface inhomogeneities, 2tt 2 Zj R 2 / L 2 ~ (_R/Aj) 2 . Therefore, as it was mentioned in 
the Introduction, the particle wavelength serves as a natural scale for describing the correlations and separates long- 
range from short-range correlations. In this sense, the interlevel transitions are more important for the short-range 
correlations. 

The diffusion coefficient is related to the mobility ( p2j ) as 

_ „ , 9 f dn^mde irti 2 L 2 h _ / R\ 

D yy = D zz = -a vy /e £ J = ^—* m = ^^<I> [z, j- j , (26) 

j 

while the mean free path along the channel is 

C = a < q > /e 2 N = (2tt) 1/2 ha TV 2 ) ^ /e 2 N 3 / 2 (27) 



IV. TRANSPORT ALONG FILMS AND CHANNELS: HIGH TEMPERATURES 

At finite temperatures all the levels with different j are populated. As a result, the transport equation becomes an 
infinite set of coupled equations. The chemical potential is the same for particles on all levels, 

where jij is the chemical potential of a 2D system of Nj fermions on level j. If we are dealing with a dilute gas, then 
[ij depends only on the number of particles on the corresponding level Zj = 2NjL 2 /ir, 

H = tfr In (l + exp (^)) = d T In (l + exp (t _ j!_ 
and 

Mi = T In (exp (jjL\ - l) (29) 



where 

2mTL 2 



d T = 



ir 2 h 2 



describes the ratio of the temperature to the energy of zero-point oscillations in the well of the width L. This equation 
should be used to express the chemical potential via the total number of particles z = 2NL 2 /ir, 
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(30) 



The solution of this equation fi (z) at T — is given by Eq. (|lfi|). 

In this Section we will calculate transport coefficients for particles with Boltzmann distribution function. In the 
case of high-temperature Boltzmann systems, 



* - #t exp (£) 9, M = T In ^ j , 6 - g exp (-2- 
and the transport equation (Q) in dimensionless variables \ji 



n j (q) = nf ) (q)[l- 



FL 3 
ir 4 T£ 2 



Xj ' 



assumes the form 



,V2 = 



1 



:{Xj(q) 



iu 2 (t/i (u, u) - ?yo (u, u)) + uj 2 (71 (u, w) - 70 (u, u)) + f (Ci (u, u) - Co (u, u)) 



2tt 2 L 2 £ 2 

S(uJ) 



(31) 



while the mobility (conductivity) is 

00 

V a® = 
Z-^i yy 

0=1 



e 2 L 2 



2tt 4 M 2 dlQ 



eE ex p(-y /" 1/2 X, (9)exp 



Here 



u = q 



L 

n'h 



S(u,j)=lnt (u + j 2 ) 



2\V 2 



u 



/2 



du 



(32) 



, u jf = u + j - j , 



and Q , 77^, and 7, are the corresponding angular Fourier harmonics of the functions 

C(q-q') =C(<Z 2 +</ 2 -2 W 'cos^), r;(q-q') =C(q-q')[l-cos^] 2 , 

7 (q - q') = C (q - q') [i - cos <p] 



(33) 



over the angle <p. In essence, the variable u = (qL/T:h) z plays the same role as the Fermi momenta zj — ^qffl L / nhj 

for degenerate systems in the previous Section. 

In the Gaussian case @, integration in ( |3l| ) leads, of course, to the same set of equations ( |l9| ) with the only 
difference that z, should be substituted by u. The situation is again non-analytic since the summation in (|3^) for 
off-diagonal transitions over j 1 should be performed up to the value S(u,j) which is not only different for each 
j, i.e., for each equation, but also depends on momentum q and exhibits step-like jumps at certain values of u = 
q 2 (L/irh) 2 . However, this non-analyticity manifests itself more noticeably in the integrands ( [32] ) rather than in 
transport coefficients themselves which are fairly smooth. Finally, the conductivity (mobility) is equal to 



°yy = 2^ G yy = ^ U I^t,— 



3 = 1 



t 5 m 2 



(34) 



U{x,y) = 



1 



x 2 Q (x) y 2 



/ Xj («) exp 



j 2 + u 2 



du 



Function n {& T ) is plotted in Fig. 4 (R/L = 0.05) and Fig. 5 (R/L = 0.5; 1). 
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FIG. 4. Function n(x), Eq.@, for R/L 
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FIG. 5. Function Il(x), Eq.fc 
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In the Boltzmann temperature range, the diffusion coefficient can be expressed via mobility as 
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drij m de To. 
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e 2 N Tr 5 h£ 2 



while the mean free path 



(35) 



(36) 



The difference between the functions II (x) in Figs. 4, 5 by several orders of magnitude is not surprising. Since x = 
■&T ~ (L/X) (A is the particle wavelength), Fig. 4 is potted in the region L ~ A. On the other hand, y = R/L — 0.05 is 
rather small meaning that R/X <C 1. As it was explained in I (and is confirmed by the present calculation), condition 
R./X <C 1 corresponds to nearly specular quantum reflection, and, therefore, to large particle mean free paths. Thus 
the large values of II (x) in Fig. 4. In Fig. 5, y = R/L ~ R/X ~ 1. This case corresponds to the most effective scattering 
of particles by surface inhomogeneities and to the smallest values of the mean free path. 
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V. SUMMARY AND DISCUSSION 



In summary, we calculated mobility and diffusion coefficients for ballistic particles in ultra-narrow channels and 
films with random rough boundaries in conditions when the motion of particles across the films is quantized. We 
obtained explicit expressions for transport coefficients via the correlation function of surface inhomogeneities. The 
most important consequence of a discrete character of the particle spectrum for the motion across the film is the non- 
analytic low-temperature dependence of the transport coefficients on the film thickness and the density of particles. 

The form of this non-analyticity strongly depends on the correlation radius of surface inhomogeneities R. In the 
case of short-range correlations of surface inhomogeneities, the low-temperature dependence of transport coefficients 
on particle density and film thickness has a pronounced saw-like structure. The saw teeth become smaller and the 
saw-like structure gradually disappears with increasing correlation radius. Finally, for long-range correlations one gets 
only not very well pronounced kinks, instead of the saw teeth, at critical values of density and/or thickness at which 
the number of occupied levels changes by one. 

Though both the amplitude an the correlation radius of surface inhomogeneities affect the particle scattering by 
the walls, the dependence of transport coefficients on the amplitude of the surface inhomogeneities £, in contrast to 
their dependence on the correlation radius i?, is quite trivial, and reduces to a multiplicative factor l/£ 2 . 

In general, the non-analytic nature of the curves is explained by the singularities in low-temperature distribution 
of degenerate fermions over a system of discrete energy levels. However, the sharp discontinuities on the saw-like 
curves for transport coefficients are caused not by the singularities in the density of state, but mostly by the interlevel 
transitions caused by the scattering on random rough walls. The beginning of the occupation of a new level leads 
to two transport effects: to the direct transport contribution of the particles from this new level, and to the opening 
of new scattering channels for particles on all already occupied levels (interlevel transitions to and from the new 
level). The first effect is proportional to the number of particles on the new level and is small. For this reason 
the non-analyticity of the transport coefficients reduces, in the absence of interlevel transitions, to a series of kinks 
corresponding to the appearance of the new levels. On the other hand, the opening of new scattering channels 
with the interlevel transitions to and from newly occupied levels affects particies on all already occupied levels thus 
increasing dramatically the total effective scattering cross-section in a step-iike manner. If one artificially freezes 
these transitions, the transport curves will exhibit kinks rather than the saw teeth. Not surprisingly, the contribution 
of interlevel transitions depends exponentially on the ratio of the particle wavelength to the correlation radius of 
the surface inhomogeneities, and decreases rapidly with increasing correlation radius of surface roughness (i.e., with 
flattening of surface inhomogeneities). 

Note, that the parameterization of transport parameters in this paper is slightly different from I. In the case of the 
mean free path it is, probably, better to use, instead of (^7j) , (j36[), the parameterization in the form I 

C~^f(R/\) 

with the minimum at R ~ A. The transformation of the results to this form is fairly obvious in both degenerate and 
Boltzmann regions. 



VI. ACKNOWLEDGMENTS. 



This work was supported by NSF grant DMR-9412769. 



[1] A.E.Meyerovich, and S.Stepaniants, Phys.Rev. B 51, 17116 (1995) 

[2] A.E.Meyerovich, and S.Stepaniants, Phys.Rev.Lett. 73, 316 (1994) 

[3] Z.Tesanovic, M.V.Jaric, and S.Maekawa, Phys.Rev.Lett. 57, 2760 (1986) 

[4] N.Trivedi, and N.W.Ashcroft, Phys.Rev. B 38, 12298 (1988) 

[5] V.B.Sandomirskii, Sov.Phys.-JETP 25, 101 (1967) [Zh.Eksp.&Teor.Fiz. 52, 158 (1968)] 



11 



VII. APPENDIX. CLASSICAL AND SEMI-CLASSICAL MOTION ACROSS THE CHANNELS 



In the classical limit, when the distance between levels with different j becomes negligible, the above results should 
coincide with the results of classical calculations in Ref. I. The transition to the classical limit corresponds either to 
thick films and/or to the states with large quantum numbers, j 3> 1, when the interlevel transitions are accompanied 
by relatively small changes of the quantum number, 1 ~ Sj <C j. The coordinate transformation (ffl) and the effective 
Hamiltonian (^J) are, obviously, the same in classical and quantum cases. The matrix elements ( ]10|) of the effective 
bulk distortion (||) in the classical limit correspond exactly to the classical matrix elements in I. However, the classical 
analog of the squares of these matrix elements, and, therefore, the collision integral ([n]) are slightly different from 
Eqs.(1.17) , (1.18) , (1.20). The reason is that in I we used not the best approximation for the following combination of 
the ^-functions: 

S' (p x - P' x ) 5 (p a -V' x ) = ~\ [S 2 (p x - Pi)]' * ~S' (p x - p' x ) (37) 

In classical calculations, quadratic expressions similar to (^) are not very well defined. More accurate expressions 
can be obtained either by working with bell-shaped functions instead of 5- functions (i.e., by introducing some small 
dissipation in the equations) , or by looking at the discrete (quantum) limit when all the combinations of the Kronecker 
symbols So- arc well defined. Careful analysis of the classical limit for the square of the matrix elements ( |l0| ) 
demonstrates that in this particular case the term in the r.h.s. of Eq.(|37|) should be written as — -^-S' (p x — p' x ) 

rather than —^kS' (p x — p' x ) and neglected in the classical limit p x 3> %/L (i.e., j ^> 1). This leads to a more 
accurate classical analog of the transition probability ([))), 

WiP > p,) = C ~£^ He - e ' )X (38) 
2 Px 6(p x -p> x ) + ^-5"(p x -p x ) 

fi(p,p') = (q-q') ' (Px1 + p' x q!), 

than Eq.(1.18). 

The corresponding change in the classical collision integral does not result in any significant changes in the expres- 
sions for the classical transport coefficients. The only improvement should be the substitution of the functions 

d sin 



a + 4 tan 4 

in the integrands for all transport coefficients by 

d sin 



a + 4/3 tan 4 (9 + 8 tan 4 (9 



where, as in I, a(u) = (5/2) t Fi (7/2, 4, -u 2 ) / r Fi (3/2, 2, -u 2 ) , and (3 (u) is similar, (3(u) = 
(3/2) } F X (5/2,3, -v?)/ l F l (3/2,2, -u 2 ). 

This change in the analytical expressions leads to more accurate results. However, numerically, the change is not 
very significant. This small numerical change is illustrated in Figs. 6, 7 for the functions fp (x) and fg (x) which 
describe the transport coefficients and the mean free path for Fermi and Boltzmann gases: 

32 e 2 L 2 R 2 N „ . , h 
-xJb (x) , x 



Tr^ U 2 JDy " (4mT) 1/2 R 
, M 4 f exp[-x 2 z 2 /cos 2 e] dz dB 

lB[X> X J iF! (3/2,2, -z 2 ) cos 2 #a + 4/3tan 4 + 8tan 4 #' [ ' 



and 
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o = — ttt^ x Jf \x) , x = V2— 



fp(x) 



~ I 1F1 (3/2, 2, -x 2 cos 2 9) cos 2 a + 4/3 tan 4 + 8 tan 4 ( 



1 



1 



G?sin0 



(40) 



Needless to say, all the prefactors remain exactly the same as in I. The variable x describes the ratio X/R in Eq.(|3£ 
and R/X in Eq.pT 




0.5 1 1.5 2 2.5 3 

FIG. 6. Function fp(x); solid line - Eq.fea), dashed line - results I 
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FIG. 7. Function /s(x); solid line - Eq. 
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